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ABSTRACT 


Use  of  body  forces  to  inject  energy  into  a  plasma  offers 
certain  advantages  over  simple  energy  addition  by  Ohmic  heating. 
To  achieve  ever-increasing  levels  of  energy  per  unit  mass  by  this 
strategy  requires  detailed  and  thorough  understanding  of  high- 
interaction  magnetohydrodynamics  (HIMHD)  through  realistic 
computer  simulation.  Such  simulation  is  possible  by  the 
existing^ — untqtre  STD  ~ Research  Corporantio^ HIMHD  codes  provided 
they  undergo  further  validation  in  the  high-interaction  regime 
through  systematic  experiments  The  present  work  has  carried  out 
a  critical  assessment  of  several  methods  for  achieving 
high-interaction,  high-magnetic  Reynolds  number  MHD  flows.  It 
Indicates  that  continuous  flow  (as  contrasted  to  pulsed  flow) 
plasmajet-driven  MHD  devices  offer  the  greatest  advantages  and 
potential  for  validating  the  STD/MHD  codes  at  high  MHD 
interaction  over  a  wide  range  of  parameters  with  the  greatest 
confidence.  It  has  led  to  the  definition  of  specific 
plasmajet-driven  experiments,  utilizing  existing  equipment,  as 
the  most  effective  way  to  carry  out  this  task.  Kc  <  v  ,  ,  ;  ,  f 


FINAL  TECHNICAL  REPORT 


STDR-84-29 


HYPERVELOCITY  PLASMAS  WITH  STRONG  MHD  INTERACTIONS 


TABLE  OF  CONTENTS 

PAGE 

ABSTRACT . ii 

TABLE  OF  CONTENTS . iii 

LIST  OF  FIGURES . iv 

1.  INTRODUCTION  .  1 

2.  OBJECTIVES  AND  STATUS  .  5 

3.  RECOMMENDATIONS  FOR  THE  VALIDATION  EXPERIMENTS  .  9 

4.  CONTRIBUTING  PROFESSIONAL  PERSONNEL  .  12 

5.  INTERACTIONS  -  PRESENTATIONS  .  13 

REFERENCES . 14 

FIGURES . 17 

APPENDIX  A  AIAA-84-0155  "Strong  MHD  Interaction" . A-1 

APPENDIX  B  Paper  "Experimental  Magnetogasdynamic  Engine 

for  Argon,  Nitrogen,  and  Air,"  1960  B-1 


LIST  OF  FIGURES 


Page 

Fig.  1.  Schematic  of  HIMHD  Device  (generator/accelerator)  17 

Fig.  2.  Schematic  of  continuous  plasma  flow  HIMHD  facility.  18 

Fig.  3.  Photograph  of  STD  Research  Corporation's 

continuous  plasma  flow  HIMHD  facility  (ca.  1975)  19 

I 

I 


FINAL  TECHNICAL  REPORT 


STDR-84-29 


HYPER VELOCITY  PLASMAS  WITH  STRONG  MHD  INTERACTIONS 


1.  INTRODUCTION 

Emerging  Air  Force  mission  needs  require  major  advances  in 
space  power  and  propulsion  technology  AFOSR  recognition  of 
these  needs  has  resulted  in  a  new  "Space  Power  and  Propulsion 
Initiative".  Leading  candidate  devices  to  meet  these  needs,  such 
as  magnetoplasmadynamic  (MPD)  thrusters  and  magnetohydrodynamic 
(MHD)  generators,  involve  the  interaction  of  partially  ionized 
conducting  gases  (plasma)  with  strong  electromagnetic  fields.  To 
minimize  weight  and  maximize  performance  of  such  MPD/MHD  devices 
requires  strong  coupling  of  the  fluid  to  the  electromagnetic 
forces,  in  other  words,  strong  electromagnetic  interaction. 

Analysis  and  computer  simulation  of  high-velocity 
partially  ionized  plasma  under  conditions  of  strong 
electromagnetic-fluid  interactions  has  received  little  attention 
to  date.  The  internationally  recognized  numerical  simulation 
capabilities  of  STD  Research  Corporation  have  been  developed  to 
include  the  capability  of  calculating  the  behavior  of  a  fluid 
which  is  strongly  coupled  to  electromagnetic  forces.  However, 
only  limited  opportunity  to  explore  strongly  coupled  cases  has 
been  possible  under  the  STD  Research  Corporation  contracts  that 
provide  numerical  simulations  support  to  the  MHD  program  of 
DOE/NASA  and  the  Navy  Pulsed  Power  MHD  generator  project.  This 
has  resulted  from  a  combination  of  both  funding  limitations  and 
the  pressing  need  to  provide  an  understanding  of  the  ongoing  MHD 
experiments  which  to  date  have  been  limited  to  modest  strength 


magnetohydrodynamic  interactions.  Thus,  no  venturesome  and  bold 
applications  of  the  numerical  simulation  tools  of  STD  Research 
have  been  made . 


Even  in  the  devices  studied  to-date,  with  modest  levels  of 
magnetohydrodynamic  interaction,  results  of  numerical  simulations 
have  uncovered  significant  new  phenomena,  which  can  have  major 
impact  on  performance  and  can  lead  to  the  invention  of 
revolutionary  techniques  and  devices.  Such  phenomena  include 
generation  of  secondary  flows,  which  in  turn  cause  major 
nonuniformities  in  the  electrode  current  and  heat  transfer  and 
boundary  layer  separation  leading  to  stall  of  the  MHD  generator, 
stagnation  pressure  asymmetries,  short-duration  high-voltage 
spikes  and  other  electromagnetic  pulse-like  effects. 
Investigation  of  design  improvements  to  minimize  such  effects  and 
maximize  performance  can  be  accomplished  with  numerical 
simulations  at  a  small  fraction  of  the  cost  of  tests  in 
appropriate  generator  test  facilities. 

The  continuing  growth  in  computer  capability  and  the 
concurrent  continuing  reduction  in  computation  costs  will 
inevitably  broaden  the  use  of  numerical  simulation  in  engineering 
design  and  development.  Already  these  developments  have  had  a 
major  impact  on  a^'rcraft  design  approaches-  In  advanced  space 
power  and  propulsion  systems,  where  hardware  and  facilities  to 
test  at  very  high  power  will  be  very  expensive,  numerical 
simulations  to  minimize  required  testing  and  to  permit  larger 
systems  and  reliable  scaling  of  small  facility  test  data  will  be 
even  more  important. 


The  credibility  of  models  used  in  numerical  simulations 
can  only  be  validated  by  comparison  with  critical  definitive 


experiments.  The  program  initiated  by  the  present  work  is 
intended  to  provide  such  a  definitive  validation  of  the  extension 
of  the  STD  numerical  simulation  capability  into  the  range  of 
strong  MHD  interactions  and  raoderate-to-high  magnetic  Reynolds 
numbers . 

Various  types  of  experiments  could  be  utilized  for  this 
validation,  including  a  steady-state  MHD  generator  operating  into 
a  space-simulated  vacuum  exhaust  or  a  high  power  MPD  arc.  In  the 
proposal  to  AFOSR  dated  April  1982  (STDR-UP-003-82)  that  led  to 
the  work  reported  here ,  STD  Research  Corporation  proposed  a 
specific  pulsed,  explosive-driven,  MHD  (XMHD)  generator 
experiment  which  offered  a,  possibly,  minimal  cost  approach  to 
obtain  the  desired  conditions  of  strong  or  very  strong 
magnetohydrodynamic  interactions.  However,  subsequent  analysis 
of  the  capabilities  of  this  device  showed  that  it  did  not  live  up 
to  expectations  for  two  reasons: 

(a)  Although  it  does  represent  a  minimal-cost  approach 
for  achieving  significant  MHD  interactions,  because  of  the  short 
times  involved,  the  instrumentation  costs  for  data  acquisition 
beyond  the  simplest  current-voltage  history  (e.g.  measurements  of 
flow  non-uniformity)  are  prohibitive. 

(b)  The  range  of  test  conditions  such  as  magnetic 
Reynolds  numbers  and  interaction  parameters  that  can  be  achieved 
at  reasonable  cost  is  not  as  wide  as  can  be  achieved  with  other 
MHD  devices  for  the  same  cost. 

In  the  course  of  the  present  work,  a  comparison  of  various 
devices  for  achieving  high  MHD  interaction  at  relatively  high 
magnetic  Reynolds  number,  including  the  explosive-driven  MHD 
(XMHD)  generator,  proved  that  the  continuous  plasma-jet  driven 


MHD  devices  offer  the  widest  range  of  interaction  parameters  at 

the  highest  magnetic  Reynolds  numbers  (i.e.  the  highest  product 

Rm  X  S  X  I  where  S  is  the  interaction  parameter  based  on 
u  p  u 

velocity  and  Ip  is  the  interaction  parameter  based  on  pressure). 

The  resulting  validated  simulation  capability  will 
stimulate  and  accelerate  inventiveness  in  the  field  of  plasma 
devices  for  Air  Force  requirements.  It  can  help  develop 
flightweight  accelerators  producing  focused  jets  with  velocities 
of  the  order  of  10  m/sec  and  particle  densities  of  the  order  of 

25  3  . . . 

10  per  m  ;  and  flightweight  power  generators  producing  either  1 
to  10  megawatts  for  hours  (10^®  to  10^^  joules  in  10^  secs)  or 

g  Q 

100  to  1000  gigawatts  for  10-100  microseconds  (10  to  10  joules 
-5  -4 

in  10  to  10  sec)  with  high  repetition  rates. 


2.  OBJECTIVES  AND  STATUS 


The  main  objective  of  the  present  work  was  to  define 
appropriate  experiments  for  the  validation  of  the  STD  Research 
MHD  computer  codes  at  high  MHD  Interactions  (of  the  order  of 
interaction  parameters  of  10  or  higher),  where  the  emerging  Air 
Force  interests  lie  (e.g.,  aerospace  mass  accelerators  or  space 
power  generator  applications). 

To  accomplish  this  objective,  progress  was  also  required 
in  understanding  the  flow  distortions  associated  with  strong  MHD 
interactions  in  hypervelocity  plasmas. 

The  present  work  has  accomplished  the  above  objectives.  A 
detailed  formulation  of  a  sound  approach  to  achieve  the 
validation  objectives  stated  above  has  been  achieved,  and  it  has 
been  presented  in  a  technical  paper  at  the  AIAA  22nd  Aerospace 
Sciences  Meeting.  This  paper,  which  is  part  of  the  present 
report  as  Appendix  A,  describes  the  technical  results  of  the  main 
tasks  of  this  work. 

The  results  and  status  of  the  present  work  can  be 
summarized  by  reviewing  the  Figures  3  and  4  of  Appendix  A. 

Figure  3  of  Appendix  A  is  a  plot  of  the  magnetic  Reynolds 
number,  Rm,  as  a  function  of  the  interaction  parameter  based  on 
velocity,  ,  for  twelve  linear  MHD  devices.  At  least  nine  of 
these  devices  have  been  built  and  tested  experimentally  (albeit 
not  all  at  their  highest  interaction  levels). 

Figure  4  of  Appendix  A  is  a  plot  of  the  magnetic  Reynolds 
number,  Rm,  as  a  function  of  the  interaction  parameter  based  on 
pressure.  Ip,  for  the  same  devices. 


5 


The  pattern  that  emerges  from  a  study  of  Figures  3  and  4 

of  Appendix  A  is  clear:  The  products  Rm  x  S  x  I  ,  Rm  x  S  and 
^  u  p  u 

Rm  X  I  ,  as  well  as  the  parameters  S  and  I  ,  all  reach  their 

p  up 

maximum  values  with  the  Demetriades-type  (References  1  -  6) 

devices  which  consist  of  an  arc-jet  generated  plasma  passing 

through  an  MHD  section  (i.e.  a  section  where  MHD  body  forces  can 

be  applied  either  to  accelerate  or  decelerate  the  plasma)  and 

exhausting  into  a  vacuum  chamber.  Operation  of  these  devices 

ranges  from  many  seconds  (if  there  is  no  special  cooling  and 

2 

large  currents,  averaging  1,000  amperes/cm  or  higher,  are 

applied)  to  continuous  (if  magnetic  fields  are  applied  to 

generate  moderate  current  densities,  of  the  order  of  100 
2 

amperes/cm  ).  The  conductivity  of  the  plasma  in  these  devices 
ranges  from  hundreds  of  mhos  per  meter  to  tens  of  thousands  of 
mhos  per  meter  depending  on  the  arc-jet  enthalpy  and  the  use  or 
non-use  of  preionizing  electrodes,  seed,  etc.  They  can  be  used 
as  generators  or  accelerators  and  the  principal  investigator  has 
experience  with  building  dozens  of  these  devices  to  operate  in 
the  5  -  100  Kw ,  50  -  200  Kw  and  100  -  1,000  Kw  power  levels. 
These  devices  are  described  in  detail  in  References  1-13. 

Progress  in  understanding  the  theory  of  these  devices 
through  the  STD/MHD  codes  has  made  it  possible  to  seek 
quantitative  answers  and  uses  for  a  multitude  of  formerly 
little-understood  effects  such  as  the  magnetoaerothermal  effect, 
electrode  shorting,  arcing  and  erosion  at  high  interaction 
parameters,  plasmajet  deflection,  etc. 

Specifically,  the  use  of  these  devices  as  the  main  source 
of  experimental  data  makes  it  possible  to  vary  or  measure  (over  a 
very  wide  range  and  pretty  much  independently),  the  following 
quantities  (and,  where  applicable,  their  distributions): 


(1)  Composition  of  working  fluid, 

(2)  Enthalpy  of  plasma  and  gas  temperature, 

(3)  Electron  temperature, 

(4)  Magnetic  field, 

(5)  Current  density  and  total  current, 

(6)  Flow  velocity, 

(7)  Mass  flow  rate, 

(8)  Ratio  of  Hall  to  Faraday  current, 

(9)  Heat  flux  to  the  walls, 

(10)  Applied  or  induced  voltage, 

(11)  Flow  static  and  stagnation  pressure, 

(12)  Degree  of  excitation  of  various  species, 

(13)  Plasma  conductivity, 

(14)  Electrode  composition,  configuration,  erosion  and 
temperature 

(15)  Reaction  forces  on  magnet  and  other  parts  of  the 
device  due  to  acceleration  or  deceleration  (thrust 
and/or  drag)  as  well  as  many  other  parameters. 

All  the  components  for  this  plasma  flow  facility  exist  in 
storage  at  STD  Research  Corporation.  It  will  be  necessary, 
however,  to  reassemble  them.  A  schematic  of  the  reassembled 
facility  is  shown  in  Figures  1  and  2.  A  photograph  of  the 

assembled  equipment  is  shown  in  Figure  3. 

Use  of  this  continuous  plasma  flow  High  Interaction  MHD 
(HIMHD)  facility  will  provide  answers  to  the  main  objective  of 
validating  the  STD/MHD  codes  by  direct  observation  and 

measurement  of  effects  predicted  by  these  codes.  Some  of  these 

effects  are  boundary  layer  stall  and  Jet  deflection  due  to  the 

magnetoaerothermal  effect,  plasma  focusing  and  compression 
effects,  etc.  These  effects  and  their  method  of  determination  or 
measurement  are  described  in  References  6,  7,  12,  14  and  15. 


In  addition,  specific  solutions  will  be  provided  to  the 
following  problems  that  still  plague  the  engineers  and 
experimenters  in  this  area  at  strong  MHD  interaction: 

(1)  Electron  emission  and  absorption  from  materials  in 
dense  plasmas  with  and  without  magnetic  fields 
(surface  effects). 

(2)  Core  and  boundary-layer  turbulence  and  related 
effects  such  as  boundary-layer  stall  and  flow 
separation  in  MHD  (interface  and  bulk  effects). 

(3)  High  interaction  effects  in  the  core  of  the  MHD 
flow-field  such  as  secondary  flows  (bulk  effects). 

(4)  Electrode  erosion,  arcing,  shorting,  aging 
enhancement  and  current  leakage  (surface  plus 
interface  effects). 

(5)  Power  take  off,  electrode  consolidation  and 
equipotent ial  exploitation  by  boundary  layer 
compensation  in  MHD  devices  (aiming  to  get  rid  of 
electrode  consolidation  circuits). 

(6)  Optimum  conductivity  and  analytic  expressions  for  the 
power  of  continuous  MHD  devices  at  any  Rm  (as  opposed 
to  pulsed  MHD  devices). 


RECOMMENDATIONS  FOR  THE  VALIDATION  EXPERIMENTS 


At  least  two  unambiguous  and  uncompromising  experiments 
can  be  carried  out  quite  easily  with  the  available  HIMHD 
equipment  to  test  the  STD  HIMHD  computer  codes.  One  experiment 
consists  of  applying  a  magnetic  field  across  the  stabilized 
arc-generated  jet  of  plasma  (either  with  simultaneously  applied 
electric  currents  or,  by  shorting  the  electrodes  of  Figure  1, 
with  an  induced  current)  and  observing  the  angle  of  deflection  of 
the  jet.  This  deflection  angle  is  due  to  the  Hall  current  and 
related  secondary  flows  and  is  known  to  exist  (References  12,  15 
-  19  and  Appendix  A)  in  HIMHD  devices.  The  other  experiment 
consists  of  an  MHD  generator  designed  to  operate  between  the 
static  pressure  furnished  at  the  exit  of  the  plasma  jet  (between 
2  and  0.1  atmospheres)  to  the  vacuum  tank  pressure  (between  0.001 
and  0.0001  atmospheres)  and  observing  the  point  at  which  boundary 
layer  over  the  anode  stalls  (e.g.  due  to  the  magnetoaerothermal 
effect)  and  the  flow  begins  to  separate.  This  stall  is  known  to 
be  caused  by  HIMHD  effects  and  has  been  observed  (References  15 
and  17  -  19). 

An  additional  argument  for  the  second  experiment  on  gener¬ 
ator  stall  is  the  following;  It  is  alright  to  talk,  on  paper, 
about  space  applications  of  MHD  generators  with  enthalpy  extrac¬ 
tion  of  25%  or  more.  It  is  well  known  that  the  higher  the 
pressure  ratio  over  which  an  MHD  generator  works,  the  higher  the 
enthalpy  extraction.  However,  the  STD/HIMHD  codes  show  that  the 
more  the  MHD  flow  is  slowed  down  and  the  faster  it  is  slowed 
down,  the  more  likely  it  is  that  it  will  separate  through 
boundary  layer  stall.  What  is  (a)  the  rate  (e.g.  with  distance 
down  the  channel)  at  which  energy  can  be  removed  from  the  flow  in 
an  MHD  generator  and  (b)  the  amount  of  energy  that  can  be  removed 
before  the  generator  stalls?  What  is  the  optimum  trade-off 


between  energy  density  and  expansion  ratio?  Clearly,  an  infi¬ 
nitely  long  MHD  channel  is  an  absurdity.  Where  do  we  cut  it  off 
for  space  applications? 

The  calculation  of  these  two  effects  in  detail  and  the 
comparison  of  the  detailed  computations  with  the  experiments 
should  furnish  the  necessary  and  sufficient  proof  required  to  use 
the  STD/MHD  codes  with  confidence  in  high  MHD  interaction 
regimes . 

The  approach  formulated  above  for  the  validation  of  the 
STD/MHD  codes  at  high  interaction  presents  the  following 
advantages  compared  to  the  original  approach  we  had  contemplated 
in  the  proposal,  STD-UP-003-83  of  April  1983,  for  the  present 
work : 

(1)  The  facility  operating  envelope  is  wider  with  respect 
to  achievable  interaction  parameters,  magnetic 
Reynolds  numbers,  and  pressure  ratios. 

(2)  More  parameters  can  be  measured  and  more  data  points 
with  the  same  time  resolution  can  be  acquired  in  a 
continuous  facility  in  the  experiment  time  available. 

(3)  Safety  precautions  are  less  demanding. 

(4)  Visual  observations  are  more  readily  available. 

(5)  Flow  chemistry  may  be  made  as  simple  or  as  complex  as 
necessary  for  valid  simulations. 

(6)  More  efficient  use  of  available  manpower. 


The  experiments  that  we  are  proposing  here  have  been  used 
at  STD  Research  Corporation  for  many  years.  In  fact,  the  basic 
techniques  have  been  described  in  Reference  1.  For  ease  of 
reference,  we  are  attaching  Reference  1  as  Appendix  B  to  this 
report . 

The  parameters  that  can  be  measured  in  these  experiments 
with  the  instrumentation  now  available  at  STD  Research 
Corporation  are  indicated  in  References  20  -  24,  where  also  the 
available  STD  continuous  plasma  flow  facility  (STD/CPFF)  and  its 
operating  envelope  are  described. 
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INTERACTIONS  -  PRESENTATIONS 


As  mentioned  above ,  the  main  results  of  the  present 
program  were  presented  at  the  22nd  AIAA  Aerospace  Sciences 
Meeting  (AIAA  Paper  84-0155  "Strong  MHD  Interaction"  by  S.  T. 
Demetriades  and  C.  D.  Maxwell)*  This  paper  created  significant 
interaction  with  other  workers  in  this  field  who  are  increasingly 
utilizing  the  STD/MHD  code  results  to  understand  and  interpret 
the  results  of  their  experiments  (e.g.,  Stanford  University  and 
AVCO  Everett  Research  Laboratories).  In  addition,  this  paper 
pointed  out  a  common  error  in  extrapolating  low  Magnetic  Reynolds 
Number  results;  it  showed  that  the  total  power  to  the  load,  Pj^  , 
varies  like  the  tanh  which  approaches  an  asymptotic  limit  for 
Rm  >  3.  Investigators  (e.g..  Reference  25)  who  use  the  low 
magnetic  Reynolds  number  approximation  to  this  relationship  (i.e. 
Pj^~Rm)  will  be  disappointed  in  their  performance  expectations  as 
they  push  velocity  and  electrical  conductivity  to  extreme  values. 
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Abstract 

Plasoa  flows  with  strong  oagnetohydrodynaalc 
(MHD)  Interaction  are  exaalned  In  nuaerlcal 
slaulatlons  of  situations  of  Interest  for  high 
perforaance  MHD  energy  conversion.  The  laportance 
of  flow  distortions  and  nonideal  processes 
associated  with  strong  MHD  Interaction  la  shown  to 
be  well  established  experlaentally.  Strategies  for 
managing  these  processes  are  exaalned  In  order  to 
extend  present  MHO  performance  levels  beyond  the 
0.6  -  0.7  HJ/kg  specific  energy  extraction 

demonstrated  to  date  In  combustion  driven 
generators.  The  role  of  the  plaaaa  boundary 
configuration  in  determining  the  flow  development 
and  the  likelihood  of  the  development  of  the 
magnetoaerotheraal  Instability  Is  Illustrated. 
Finally,  the  Importance  of  high  ugnetlc  Reynolds 
number  effects  are  considered  for  MHD  devices 
capable  of  high  levels  of  MHD  Interaction. 

Introduction 

We  consider  here  effects  Inherent  la 
high  performance  aagnetohydrodynaalc  energy 
conversion  devices,  particularly  those  which  employ 
plasma  as  the  worUng  fluid.  Five  nondlaanalonal 
parameters  characterize  the  global  aagnetohydro- 
dynamics  in  such  devices;  the  Mach  number,  an  MHD 
Interaction  parameter  defined  below,  the  aspect 
ratio  (L/D),  the  viscous  Reynolds  number.  Re,  and 
the  magnetic  Reynolds  number,  Rm.  Ihimerous  local 
magnetogasdynaalc ,  thermal,  kinetic,  and  electrical 
nondlmenslonal  parameters  also  govern  the  local 
behavior  of  such  devices. 

Given  practical  constraints  on  magnetic  field 
(1-10  T),  length-co-dlameter  ratio  (10  to.lOlX}), 
and  mass  density  of  the  plasma  (  >  0.01  kg/m^).  such 
devices  may  have  length  scales  of  the  order  of  0.1 
to  100  m.  Figs.  1  and  2  give  typical  values  of 
these  constraints  and  the  feasible  domain  for  two 
MHD  accelerator  devices .  The  Mach  ouad>er  will 
vary  from  high  subsonic  for  some  fossil-fueled 
devices  to  M  >  10  for  some  explosive  devices.  ,The 
viscous  Reynolds  number  will  be  of  the  order  10^. 

While  Idealized ,  one-dimensional  models  of  the 
flow  In  such  devices  may  be  adequate  to  describe 
situations  In  which  the  Interaction  between  the 
working  fluid  and  the  ponderomotlve  (Lorentz) 
forces  Is  weak,  more  elaborate  models  must  be 
employed  to  simulate  devices  In  which  the 
Interaction  parameters  and/or  magnetic  Reynolds 
number  reach  high  values.  Highly  unconventional 
and  three-dimensional  flow  behavior  results  from 
strong  Interaction  between  the  Lorentz  forces  and 
the  fluid.  This  paper  reviews  the  operating 
regimes  of  major  classes  of  MHD  power  generators 
tested  to  date ,  examines  unconventional  phenomena 
predicted  to  occur  In  such  devices,  and  describes 
the  recent  experimental  confirmation  of  the 
magnetoaerothermal  effects  predicted  by  the  STD/MHD 
codes  and  reported  for  the  first  time  three  years 
ago  at  this  meeting. 

*Presldcnt  and  Technical  Director,  Associate  Fellow  AIAA 
tvice  President  and  MHD  Project  Manager,  Menber  AIAA 


The  STD/MHD  Codes 

The  physical  processes  which  govern  MHD 
device  performance  are  numerous  and  Interdependent , 
and  crucial  phenomena  nay  easily  be  obscured  by 
intuitive,  simplified,  or  decoupled  analyses.  Many 
such  processes  are  Inherently  unsteady,  probabil¬ 
istic,  and/or  three  dimensional. 

Comprehensive  analytical  tools  are  required  for 
unambiguous  Interpretation  of  MHO  experimental  data 
and  to  extend  test  .  data  to  new  designs,  larger 
scales,  or  higher  levels  of  MHD  Interaction.  The 
STD/MHD  codes  are  constructed  from  a  large 
collection  of  modules  or  subprograms  which  address 
various  aspects  of  the  MHD  problem.  Taken  to¬ 
gether,  these  codes  comprehensively  define  and 
solve  the  MHD  problem. 

The  theoretical  foundation  for  these  modules 
has  been  laid  over  the  past  25  years  by  the  senior 
author  and  his  coworkers.  The  formulations  for 
most  of  these  modules  are  documented  In  Ref .  [11 
and  Its  references,  as  well  as  several  more  recent 
STD  Research  Corporation  publications  (e.g. ,  [2]  - 

C*l). 

Post-test  simulations  of  well-diagnosed 
experiments  have  never  failed  to  verify  the 
accuracy  of  the  STD/MHD  codes  within  the  experi¬ 
mental  tolerances  of  the  measured  Input  data  [5], 
[6].  The  STD/MHD  codes  have  been  shown  to  be 
sufficiently  sensitive  to  detect  subtle  Inconsis¬ 
tencies  In  Input  data  furnished  by  experimenters. 
These  Inconsl.  tencles  were  confirmed  by  the 
experimenters  only  after  they  were  pointed  out  to 
them  by  STD  Research  [5].  The  STD/MHD  codes  have 
been  used  as  engineering  tools  to  design  MUD 
devices.  Improve  the  performance  of  devices 
designed  by  others,  predict  nominal  performance, 
and  assess  the  Impact  of  design  changes  for  a 
number  of  new  MHD  devices .  These  designs ,  Improve¬ 
ments  and  predictions  have  always  been  proven  or 
verified  when  the  devices  were  tested  [6]. 

Host  importantly.  In  a  number  of  cases 
anomalous  or  unexpected  phenomena  have  emerged  from 
three-dimensional  and/or  time-dependent  slanilatlons 
with  the  STD/MHD  codes  before  experimental 
observations  of  the  phenomena  were  available. 
Many  of  these  effects  have  been  verified 
subsequently,  and  none  have  been  dlsproven.  For 
many  processes,  such  as  the  magnetoaerothermal 
effect  [7|,  the  phenomena  Involve  such  complex 
Interactions  that  It  Is  doubtful  that  a  prediction 
or  complete  understanding  would  have  been  possible 
without  the  physical  Insight  provided  by 
comprehensive  numerical  simulations  and  adequate 
computer  graphics. 

These  codes  are  available  and  accessible  to  all 
who  need  to  use  them. 
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MHD  PerfortBance  and  MHD  Interaction 

While  the  STD/MHD  codes  have  also  been 
validated  in  many  weak-interaction  MHO  situations, 
(indeed,  with  no  magnetic  field  they  have  been 
validated  as  computational  fluid  mechanics  codes), 
it  is  at  moderate  to  high  levels  of  MHO  interaction 
that  a  comprehensive,  fully  coupled  analytical 
approach  is  required.  Since  it  la  also  at  high 
levels  of  MHO  interaction  that  the  performance  of 
MHD  devices  becomes  attractive,  simplified 
treatments  have  limited  usefulness,  even  for  system 
calculations. 

The  relationship  between  performance  and  MHO 
interaction  is  well  known.  One  way  of  expressing 
this  relationship  [8]  is  that  the  energy-per- 
unlt-mass,  or  specific  energy,  that  la  added  or 
removed  from  the  plasma  flow  in  MHD  devices  varies 
directly  as  the  product  of  the  conversion 
efficiency  and  Che  interaction  parameter  I  , 


where  we  have  used  Che  symbol  P  for  the  electrical 
power  and  m  for  Che  mass  flow  rate.  The  propor¬ 
tionality  constant  in  this  expression  is  approxi¬ 
mately  RI.  the  product  of  Che  gas  constant  and  the 
temperature  of  the  working  fluid.  The  conve»^on 
efficiency  is  defined  by  n  *  7-T  /  7^(t)xB). 
while  Che  interaction  parameter  I  is  defined  over 
the  length  L  of  an  MHO  device  by  tfte  expression 

I  .  /  llill  dx  (2) 

OP 

where  p  is  the  static  pressure  in  the  flow. 

It  is  noted  chat  ocher  equivalent  forms  of  Eq. 
(1)  may  be  written.  These  may  Involve  other 
measures  of  performance,  such  as  enthalpy 
extraction  ratio.  However,  it  is  not  strictly 
valid  CO  write,  as  some  have  done,  performance  as 
varying  with  the  interaction  parameter  I  alone;  to 
cite  an  extreme  example,  high  values  of"  I  may  be 
reached  by  short-circuiting  a  generator,^  but  no 
power  is  produced.  Performance  requires 
interaction  and  efficiency. 

One  can  also  consider  two  Interaction 
parameters,  namely,  the  ratio  of  Lorentz  force  per 
unit  volume  to  the  mom^cum  flux  or  dynamic 
pressure,  Ij  ■  pU^),  and  the  ratio  of 

Lorentz  force  per  unn  volume  to  static  pressure  p, 
i^  *  (J  B/p).  These,  as  well  as  the  local  magnetic 
Reynolds  number,  r  •  (rU,  can  be  made 

nondlmensional  by  integrating  over  a  length,  L. 
The  nondlmensional  interaction  parameter  obtained 
by  integrating  i^  is  called  Che  interaction 
parameter  based  on  velocity,  S  .  The 

nondlmensional  interaction  parameter  (Stained  by 
integrating  1^  is  the  Interaction  parameter  based 
on  pressure,  I  ,  defined  above.  Rm  is  Che 
nondlmensional  m^necic  Reynolds  number. 

2 

It  is  noted  that  1,  ■  l,(c  /c^)M  ,  where  c  /c 
is  Che  ratio  of  specific  heats  and  M  is  Che  RacK 
number.  Thus  for  a  given  performance  level 
characterized  by  I  ,  S  will  vary  inversely  as  the 
square  of  the  chjracc^erisclc  Mach  number  in  the 
device.  While  1  correlates  performance,  S 


correlates  flow  distortion  and  behavior  (8|.  It 
follows  chat  for  two  devices  with  Che  same 
performance,  expressed  in  power  per  unit  mass  flow, 
Che  device  with  the  lower  Mach  number  will  exhibit 
Che  greater  flow  distortion. 

Figs.  3  and  4  present  a  map  of  Che  character¬ 
istic  parameters  I  ,  S  and  Rm  for  MHO  devices 
using  as  working  !lul(^  combustion  products  of 
fossil  fuels  and  solid  propellants,  explosion 
products  and  noble  gases,  and  arc-jet  heated 
nitrogen  and  argon.  The  highest  levels  of 
interaction  were  achieved  in  the  arc-jet  experi¬ 
ments,  which  were  run  both  in  the  generator  and 
accelerator  modes  t91  -  [Uj.  The  highest 
magnetic  Reynolds  numbers  were  obtained  in  the 
explosive  MHO  experiments,  %ihlch  also  exhibit  the 
lowest  values  of  because  of  high  Mach  numbers. 

The  Appendix  Co  this  paper ,  taken  from  [  12  ] , 
derives  expressions  for  the  power  from  an  idealized 
MHD  generator  with  arbitrary  magnetic  Reynolds 
number.  It  is  shown  Chat  power  varies  as 
tanh(Rm/2).  For  Rm  «  1,  the  power  varies  as  Rm  in 
accordance  with  textbook  formulas  for  low  magnetic 
Reynolds  number  HUD.  Ho«fever ,  for  Rm  of  order 
unity  or  greater,  power  tends  coward  an  asymptotic 
limit  as  Che  hyperbolic  tangent  tends  coward  1. 
For  Rm  •  3,  power  has  reached  90Z  of  its  asymptotic 
value.  Further  Increases  in  conductivity  and 
velocity  for  the  sake  of  increasing  power  yield 
diminishing  returns,  especially  when  losses  are 
considered. 

The  validity  and  universality  of  the  correl¬ 
ation  in  Eq.  (1)  is  shown  by  Fig.  5.  The  specific 
energy  extraction  P/m  is  plotted  versus  the  product 
T)  I  for  fossil-fueled  and  solid  propellant  experi¬ 
ments  and  designs  in  the  U.S.  and  elsewhere.  These 
points  have  all  been  plotted  tollowlng 
determination  of  these  parameters  from  the 
experimental  data  or  simulations  of  proposed 
devices  with  the  STD/MHO  codes.  These  codes  auto¬ 
matically  display  such  parameters  as  n  and  1  as 
well  as  the  extracted  electrical  power.  ^The 
experiments  span  a  wide  range  of  length  scales  (O.S 
CO  30  m) ,  configurations  (linear  or  disk,  square, 
octagonal,  round,  etc.),  electrical  connections 
(Hall,  Faraday,  wlndowframe  and  external  diagonal), 
and  working  fluids  (toluene,  JP-4,  coal,  natural 
gas,  solid  propellants,  etc.,  all  at  various 
stoichiometries ) . 

The  levels  of  interaction  spanned  in  Fig.  5 
have  been  subdivided  into  four  regimes  based  on  a 
qualitative  assessment  of  the  degree  of  MHO  flow 
distortion  observed  in  the  simulations  of  these 
devices.  Because  the  degree  of  MHO  flow  distortion 
at  a  given  I  also  depends  upon  the  Mach  number  in 
each  device ,  ^these  suMivlslons  are  approximate. 

Many  experiments  have  been  carried  out  in  the 
low-interaction  regime,  where  the  fluid  behavior  is 
conventional,  i.e.  as  in  ordinary  fluid  mechanics, 
with  only  second-order  corrections  for  MHD  effects. 
The  simulation  of  these  experiments  with  Che 
STD/MHD  codes  has  shown  chat  Che  second  order 
modifications  Co  friction  and  heat  transfer  due  to 
the  weak  Lorentz  forces  and  ohmic  heating  alter  the 
conventional  duct-flow  behavior  very  slightly.  The 
weak  secondary  flow  pattern  arising  from  the 
viscous  interaction  between  the  flow  and  Che  walls 
(an  eight  cell  pattern  in  a  square  duct)  may  be 
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disrupted  by  Che  weak  MHD  InCeracclona ,  which 
proance  nonunlfom  transverse  body  forces,  but  the 
bulk  flow  properties  are  virtually  unaffected.  The 
flow  field  can  then  be  adequately  predicted  with 
the  isethods  of  conventional  gasdynamlcs  aerely  by 
accounting  for  Che  additional  average  axial  body 
forces  and  Che  bulk  oholc  heating  of  the  gas. 

The  few  experiments  that  have  been  carried  out 
In  the  high-lnteracclon  regime  (see  Fig.  S)  have 
shown  significant  MHD  effects  on  the  bulk  flow  of 
the  plasma,  notably  flow  asymmetry  and  distortion. 
These  hlgh-lnteractlon  experiments  Include  the  AVCO 
Mark  V  (13],  the  Viking  I  teats  at  the  Wright 
Patterson  APL  (14),  the  Bangerter  solid  propellant 
experiment  at  Hercules,  Inc.  [IS],  the  Pamir  I 
experiment  with  solid  propellant  (16),  and  Che  AEOC 
High  Performance  Demonstration  Experiment  (HPDE) 
(17).  In  all  of  these  experiments  the  specific 
energy  extraction  reached  at  least  0.6  MJ/kg.  STD 
Research  Corporation  has  performed  extensive 
simulations  of  these  experlawncs,  especially  the 
AEDC/HPDE,  and  has  shown  the  Influence  of  these 
effects  In  detail ,  as  will  be  discussed  below. 

The  importance  of  hlgh-lnteractlon  phenomena 
has  been  appreciated  by  many  of  Che  experiment¬ 
alists  who  have  achieved  the  highest  levels  of 
performance .  For  example ,  at  the  conclusion  of  the 
Pamir  I  experiments  In  the  Soviet  Union  (16)  it  was 
noted  chat , 

'Using  an  effective  solid-fuel  plasma 
generator,  specific  power  output  and  power 
density  of  up  to  0.6  MJ/kg  and  500  W/cm  , 
respectively,,  was  generated.  Further 
Increasing  these  parameters  (and  accordingly 
...  the  coefficient  of  enthalpy  extraction)  can 
be  achieved  only  If  the  limiting  effect  of 
processes  associated  with  strong  interaction 
are  suppressed.  The  nature  of  these  processes 
Is  not  yet  fully  understood,  thus  It  Is 
necessary  to  continue  Che  study  of  these 
phenomena . ' 

In  Che  high  and  very  high  Interaction  regimes, 
where  the  nonlinear  coupling  of  Che  flow  field  to 
the  Lorentz  forces  and  the  ohmic  heating  will  be 
very  strong,  the  plasma  flow  will  be  dominated  by 
MHD  effects.  This  Is  clearly  Indicated  by  the 
experience  obtained  by  the  STD/MHD  code  simulations 
as  the  Interaction  level  and  performance  Increase. 
This  experience  will  be  reviewed  In  the  following 
section. 


High  Interaction  MHD  Phenomena 

The  successful  validation  and  use  of  the 
STD/MHD  codes  under  moderate-co-hlgh  levels  of  MHD 
Interaction  have  demonstrated  chat,  as  the  MHD 
Interaction  Increases,  profound  changes  can  cake 
place  In  the  fluid  dynamic  behavior  of  the  duct. 
The  following  unconventional  effects  are  some  of 
chose  that  have  been  observed  In  both  the  computer 
simulations  and  the  experimental  data  obtained  to 
date: 

(1)  Sidewall  velocity  overshoots.  Diminished 
axial  Lorentz  forces  In  the  boundary  layers  on 
Insulating  sidewalls  permit  the  sidewall 
working  fluid  to  accelerate  relative  to  the 
core  fluid  [18],  [8].  When  current  Is  allowed 


to  flow  to  Che  sidewalls  of  a  generator,  the 
resulting  decrease  In  Lorentz  force  over  the 
electrodes  can  lead  to  similar  velocity 
overshoots  In  the  electrode  boundary  layers 
[19].  In  channels  of  high  Interaction,  neglect 
of  this  effect  leads  to  serious  errors  In 
determining  the  blockage  at  the  diffuser 
entrance  and  the  resulting  diffuser  performance 
[81. 

(2)  Transverse  static  pressure  anoMlles.  The 
virtually  unavoidable  presence  of  local  axial 
currents  promces  spatially  nonuniform  body 
forces  on  Che  working  fluid  In  Che  direction 
perpendicular  to  both  the  mean  flow  and 
magnetic  field  directions.  In  the  case  In 
which  the  average  of  these  transverse  body 
forces  Is  not  zero  In  the  cross  section,  a 
transverse  pressure  gradient  Is  sustained  [19]. 
This  nay  occur  even  In  a  segmented  Faraday 
device  In  which  the  net  Hall  current  Is  zero. 

(3)  Transverse  asymmetries  In  the  electric 
field  distribution.  Secondary  flows  and  the 
magnetoaerothermal  effect  redistribute  the 
momentum  and  energy  fields  of  the  working 
fluid.  This  Is  observed  in  measurements  of  the 
transverse  voltage  profiles  of  moderace-co-hlgh 
Interaction  MHD  generators  [6]. 

(4)  Transverse  asymmetries  In  the  velocity 
field.  Nonuniform  Lorentz  force  distribution 
In  the  cross  section  results  in  skewed  profiles 
of  axial  velocity.  Such  velocity  asyasetrles 
are  manifested  In  asymmetric  stagnation 
pressure  profile  measurements  [21]. 

(5)  The  magnetoaerothermal  effect.  This 
effect  causes  asymmetric  boundary  layer 
deceleration  or  stall  (separation)  at  the 
center  of  the  anodes  In  an  MHD  generator. 
Even  In  segmented  Faraday  generators  or 
accelerators  with  significant  Interaction 
levels,  correct  analysis  shotw  that  there  exist 
significant  Hall  currents  nonunlformly 
distributed  In  the  cross  section,  even  at  zero 
Integrated  Hall  current.  The  nonuniformity  In 
the  transverse  Lorentz  force  field  caused  by 
these  Hall  currents  drives  strong  secondary 
flows  (In  confined  channel  flom  what  goes  up 
must  come  down;  In  free  Jets  the  Jet  Is 
permanently  deflected;  under  the  right 
conditions  the  flow  oscillates)  [7]  and  [22]. 
The  secondary  flow  components  on  the  transverse 
plane  direct  hot  plasma  on  the  center  of  the 
anode  walls.  This,  among  other  things, 
constricts  the  discharge  or  current  paths  over 
the  anode  (which  promotes  stronger  Hall  current 
gradients  and  secondary  flows),  causes  a  large 
Increase  In  the  local  current  density,  and 
Increases  the  local  Interaction  (hence  the 
local  or  not-so-local  stall  In  generators)  and 
the  local  heat  transfer.  The  cathode  boundary 
layer  Is  also  Influenced  by  the 
magnetoaerothermal  effect  and  may  show  a 
tendency  to  stall  further  downstream  and  away 
from  the  cathode  center  In  segmented  Faraday 
generators . 


Fig.  4  illustrates  these  effects,  especially 
(1),  (3),  (4),  and  (5),  on  the  basis  of  the  STD/MHD 
code  simulations  of  the  AEDC/HPDE  Run  MI-006-014. 
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This  figure  gives  a  capsule  view  of  the  effects  of 
HHD  Interaction.  Because  the  Interaction  paraaeter 
Increases  directly  with  length,  the  velocity  and 
teaperature  profiles  and  current  and  secondary  flow 
patterns  at  successive  axial  stations  Illustrate 
the  aagnltude  of  HHD  flow  distortion  found  at 
successively  higher  levels  of  HHD  InCeractioa. 

To  Illustrate  how  the  aagnetlc  field  and  wall 
teaperature  can  Influence  the  development  of  high 
Interaction  HHD  phenoaena.  Figs.  7  and  8  show  a 
progression  of  computations  with  peak  aagnetlc 
field  Increasing  from  3  to  S  tesla  for  the 
conditions  of  the  AEDC/HPDE  Bun  Ml-006-014  (see 
Fig.  4),  except  that  the  electrode  wall  teaperature 
has  been  elevated  from  values  of  1050  K  aaxlaua  to 
values  of  1700  K  aaxlaua.  The  developaant  of  the 
velocity  profiles  and  secondary  flows  shown  In  Fig. 
7  and  the  developaant  of  the  static  teaperature  and 
transverse  current  distribution  shown  In  Fig.  8 
clearly  deaonstrate  the  strengthening  of  the  above 
effects  as  aagnetlc  field  Increases.  Equally 
important  Is  the  observation  that  the  strength  of 
the  above  effects  Is  dlalnlshed  by  the  Increase  In 
wall  teaperature,  as  was  shown  In  Ref.  (22|. 
Hotter  walls  proaote  a  more  lailfora  conductivity 
field,  which  leads  to  weaker  Hall  current  gradients 
In  the  cross  section  [7|. 

All  of  the  above  effects  are  now  well 
established  experlaentally : 

(1)  Sidewall  velocity  overshoots  have  been 
directly  aeasured  by  Rankin  at  Stanford 
University  (23|  and  more  recently  by  McClalne, 
et  al.  at  Avco  Everett  Research  Laboratory 
T?41. 

(2)  Transverse  static  pressure  gradients  have 
been  observed  In  every  high  interaction 
experlaent  In  which  they  have  been 
Investigated.  Especially  severe  static 
pressure  anomalies  are  observed  In  channels 
which  permit  a  significant  net  Hall  currant  to 
flow  (e.g.,  ( 16|). 

(3)  Review  of  the  transverse  voltage  profile 
measurements  over  the  past  20  years  reveals  a 
qualitatively  similar  asyametry  In  them  which 
la  due  at  least  In  part  to  the  redistribution 
of  momentum  and  energy  In  the  cross  section  by 
HHD  Interactions.  This  asymmetric  shape  Is  a 
characteristic  of  all  of  the  hlgh-lnteraction 
three  dimensional  simulations  performed  to 
date ,  and  comparisons  between  three  dimensional 
simulations  and  experimental  profiles  have 
yielded  very  satisfactory  agreeme  nt  (6|. 

(4)  Velocity  profiles  Inferred  from  exit 
stagnation  pressure  profile  measurements  have 
been  shown  to  be  Increasingly  unconventional 
and  asynsetrlc  from  anode  to  cathode  as  ^BlD 
Interaction  Increases  [21],  [24].  While  there 
has  been  a  cautious  tendency  by  the 
experimenters  to  look  for  possible  explanations 
for  these  asymmetries  In  conventional  effects 
(e.g.,  shocks)  or  errors  in  experimental 
technique  (e.g.,  probe  Interactions)  [24],  the 
trends  in  these  measurements  are  now  believed 
to  be  sufficiently  systematic  and  consistent 
with  predictions  of  the  analytical  models  to 
conclude  that  these  asymmetries  represent  high 
interaction  HHD  phenomena. 


(5)  The  above  effects  (2)  -  (S)  are  closely 
related  to  Che  magnecoaerochermal  effect.  In 
fact ,  they  may  be  thought  of  as  precursors  or 
signatures  of  an  Incipient  magnecoaerochermal 
current  conscrlctlon/anode  boundary  layer 
stall.  Following  1981  predictions  by  STD 
Research  Corporation  chat  the  magnecoaero¬ 
chermal  effect  should  be  observed  In  the 
AEDC/HPDE  tests,  such  a  magneocaerochermal 
effect  was  clearly  observed  In  Run  HI-007-016 
of  20  October  1982.  The  follo%>lng  section 
describes  the  AEDC/HPDE  simulations  and  data 
related  to  the  prediction  and  subsequent 
observation  of  the  magnecoaerochermal  effect. 


The  most  recent  and  the  most  convincing 
confirmation  of  the  Importance  of  the  above  effects 
has  come  In  the  careful  and  well-diagnosed 
AEDC/HPDE  tests  carried  out  over  the  past  three 
years. 

Prior  to  Che  first  runs  of  the  HPDE,  the  three 
dimensional  simulations  carried  out  for  a  variety 
of  channels  from  low  to  high  Interaction  showed  the 
tendency  of  the  anode  boundary  layer  to  be  stalled 
by  a  current  concentration  In  Its  center  at  high 
Interaction.  Pretest  simulations  at  the  nominal 
design  conditions  of  the  HPDE  clearly  showed  the 
presence  of  this  behavior.  Depending  on  wall  and 
operating  conditions,  the  effects  also  appeared  In 
simulations  with  less  than  the  nominal  aagnetlc 
field  strength  of  6  tesla.  These  results  were 
reported  In  Refs.  [7]  and  [22]. 

On  Che  basis  of  the  three-dimensional  slnula- 
Clons  of  the  AEDC/HPDE  carried  out  by  STD  Research 
Corporation,  it  was  possible  to  predict  in  early 
1981  that  [7], 


’Effects  of  Che  kind  described  above  [i.e.,  the 
aagnetoaerochemal  effects)  should  be  visible 


in  some  of  the  machines  currently  under  test . 
For  example,  the  AEDC/HPDE  should  exhibit 
accelerated  breakdown  and  erosion  of  the  center 
of  the  anodes  In  the  downstream  half  of  the 
machine .  .  . ’ 


The  first  sign  chat  the  predicted  high 
Interaction  phenomena  were  present  In  the  HPDE 
tests  was  the  measurement  of  unmistakable 
transverse  voltage  asymmetries  In  the  HI-006  test 
series  data  (1979-1980).  While  the  resolution  of 
these  data  was  then  somewhat  coarse,  the  agreement 
with  the  predictions  of  the  STD/HHD  codes  was  quite 
satisfactory  [6],  The  higher  resolution 
transverse  profiles  of  the  later  HI-007  test  series 
(1981-1982)  amply  confirmed  the  presence  of  the 
voltage  asymmetries  [4], 


There  is  a  striking  similarity  between  the 
asymmetric  voltage  profiles  in  the  AEDC/HPDE  [6] 
and  the  profiles  at  the  downstream  electrodes  in 
previous  hlgh-lnteraction  Faraday  generator  tests 
(e.g.,  the  Viking  I  tests  —  see  Fig.  32  of  Ref. 
[14]).  Some  of  the  data  from  low-interaction 
Faraday  tests  appear  to  contradict  these  trends 
(25).  However,  (a)  the  observed  differences 
between  anode  and  cathode  are  much  smaller  In  the 
low-interaction  testa  than  In  the  hlgh-lnteraction 
tests,  and  (b)  the  weak  Interaction  data  of  [25] 
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was  taken  In  the  presence  of  slag,  introducing  net 
Hall  currents  In  the  plasma  due  to  current  leakage 
in  the  slag  as  well  as  distortions  of  the  cathode 
voltages  by  the  "stepping"  or  "staircase"  slag 
shorting  effect  [26],  The  pressure  trace  of 
magnetoaerothermal  effects  is  clear  in  low- 
interaction  machines,  although  the  voltage  drop 
asyoaetry  is  not  as  pronounced.  One  thing  can  be 
said  with  certainty:  one  needs  large  scale  MHS 

channels  or  artificially  enhanced  HHD  interactions 
(e.g.,  by  artificially  increasing  the  conductivity 
and/or  the  magnetic  field  or  dropping  the  pressure, 
etc.)  to  make  meaningful,  large  slgnal-to-noise , 
measurements  of  the  magnetoaerothermal  effects. 

Beginning  with  the  Ml-007  series  of  tests  at 
the  HPDE,  the  stagnation  pressure  profile  data 
approximately  one  meter  downstream  of  the  channel 
exit  strongly  suggested  the  increasing  presence  of 
stall  on  the  anode  wall  as  the  magnetic  field  was 
increased  in  successive  runs.  Fig.  9  shows  Che 
simulated  development  of  the  stagnation  pressure 
profiles  down  the  channel  under  the  S  tesla 
conditions  of  Figs.  7  and  8.  In  Fig.  9  Che  3  m  and 
4  m  stations  have  been  expanded  and  cut  away  to 
reveal  the  centerline  total  pressure  profiles  i^ich 
were  measured  in  the  MI-007  tests.  The  (3  and  4  m 
stations  have  approximately  the  same  interaction 
parameter  as  the  HI-007  tests  due  to  Che 
above-noted  dependence  of  interaction  parameter  on 
length. )  The  qualitative  similarity  between  the 
computed  profiles  and  the  measured  profiles  Is 
quite  good. 

Run  HI-007-016  on  20  October  1982  was  carried 
out  with  Che  highest  magnetic  field  in  the  test 
series,  3.8  tesla.  For  the  first  half  of  this  7 
second  run,  Che  electrical  performance  of  the 
downstream  half  of  the  channel  wss  anomalously 
degraded.  Post-teat  Inspection  of  the  channel 
showed  significant  erosion  and  discoloration  in  a 
long  stripe  down  the  center  of  the  anode  wall  in 
the  downstream  half  of  Che  channel,  as  predicted  in 
1981.  Fig.  10  presents  a  photograph  of  this 
manifestation  of  the  magnetoaerothermal  effect. 

Significantly,  the  anomalous  electrical 
behavior  in  Run  HI-007-016  "healed"  approximately 
half  way  through  the  run.  It  is  believed  chat  this 
behavior  was  connected  with  the  rise  in  the  wall 
temperature  during  Che  test. 

Control  of  Magnetoaerothermal  Effects 

At  high  levels  of  HHD  Interaction  magneCoaero- 
chermal  effects  promote  nonunif ormicy  in  the 
electrical  and  thermal  stresses  on  the  walls. 
Therefore,  in  planning  duration  tests  under 
"cosmercial  conditions",  the  local  conditions  on 
Che  perimeter  must  be  matched,  rather  chan 
"average"  or  "core"  conditions.  The  role  of 
magnetoaerothermal  effects  in  determining  ultimate 
performance  at  the  very  high  HHD  interaction  levels 
indicated  by  simplified  system  analyses  remains  to 
be  investigated  experimentally. 

Strategies  which  can  control  and  even  exploit 
Che  magnetoaerothermal  effects  have  been 
investigated  by  STD  Research  Corporation,  [71,[22|: 

(1)  Hotter  wall  surface  temperatures  increase 
plasma  uniformity  and  diminish  the  local  Hall 
currents  which  drive  magnetoaerothermal  effects. 


(2)  For  channels  with  the  same  performance, 
higher  Mach  number  designs  minimize  the  interaction 
parameter  based  on  velocity  and  therefore  the 
strength  of  the  magnetoaerothermal  effects. 

(3)  Controlled  application  of  axial  currents 
though  the  choice  of  the  electrical  connections  and 
load  network  can  allow  the  designer  to  "steer" 
secondary  flows  and  drive  magnetoaerothermal 
effects  in  beneficial  ways. 

(4)  Application  of  inlet  swirl  can  have  a 
strong  Influence  on  the  development  of 
magnetoaerothermal  effects. 

(5)  Alternate  electrode  configurations  and 
cross  sectional  geometries  can  modify  the 
development  of  body  forces  which  drive 
magnetoaerothermal  effects.  Current  control 
techniques  whereby  current  is  forced  to  follow 
planned  paths  to  the  electrodes  slso  fall  in  this 
category. 

The  success  of  the  above  strategies  will  make 
or  break  the  design  of  high  performance  HHD  power 
generators.  We  have  been  searching  for  improve¬ 
ments  in  the  effective  uniformity  through  use  of 
nonrectangular  MHO  channel  cross  sections. 
Preliminary  results  for  a  1000  MU  (thermal)  design 
indicate  that:  (1)  optimization  of  the  channel 
shape  and  the  extent  of  electrode  coverage  on  the 
channel  perimeter  can  yield  performance  within  a 
few  percent  of  the  ideal  (perfectly  uniform) 
conditions  (Figs.  U  and  12),  and  (2)  even  in 
channels  with  circular  cross  sections,  substantial 
secondary  flows  are  generated  by  the  transverse  MUD 
body  forces  in  the  plasma  (Figs.  13  and  14);  thus 
it  Is  likely  magnetoaerothermal  effects  will  have 
to  be  reckoned  with  regardless  of  channel  shape. 

Conclusions 

This  psper  presents  detailed  results  of  the 
application  of  the  STD/MHD  codes  to  the  design  or 
analysis  of  a  number  of  specific  HHD  power 
generator  devices  chat  have  been  built  or  examined 
over  Che  years.  In  addition,  it  presents  the 
detailed  description  of  a  number  of  more  complex 
(but  equally  important)  HHD  phenomena  revealed  in 
detail,  again  for  the  first  time,  by  the 
application  of  the  STD/MHD  codes  to  specific 
devices  capable  of  operating  at  HHD  interaction  and 
magnetic  Reynolds  number  levels  approaching  the 
levels  expected  to  be  encountered  in  real 
(practical,  useful,  commercial,  etc.)  HHD 
generators. 

The  STD/MHD  codes  have  been  validated  both  by 
direct  simulation  of  carefully  diagnosed 
experiments  and,  even  more  importantly,  by  their 
demonstrated  capability  to  predict  and  interpret 
previously  unsuspected  HHD  phenomena  which  are 
crucially  important  for  the  successful  operation  of 
HHD  generators.  The  available  experimental  data 
range  from  the  low  to  the  moderate-to-high 
Interaction  regimes.  While  the  codes  have  also 
been  exercised  in  the  high-to-very  high  regimes 
desirable  for  full-scale  commercial  and  other  high 
power  systems,  there  is  a  scarcity  of  detailed  and 
sufficient  data  with  which  to  compare  the 
predictions  of  the  codes  at  such  levels  of  HHD 
interaction. 
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fcpp«ndix 

InfXu«nc«  of  Magnetic  R«voold»  MuabT  on 
Pow«r  g«nTaf  d  by  an  Id«*X  WHO  Davic« 

In  tius  Appendix  w«  dative  a  siapla  expraeaton 
Cor  the  power  oC  «n  rdeaX,  one'diaansioneX  HHO  9en- 
eracor  that  is  not  restricted  by  the  usual  assuap- 
txon  of  a  vary  nail  aagnatic  Reynolds  niatber*  It 
applies  to  pulsed  MHO  generators  and,  possibly,  to 
other  devices.  It  shows  that  as  the  aagnatic 
Reynolds  nuabar  increases,  the  power  to  the  load 
approaches  an  asyaptotic  liait, 

Consider  the  MHO  problaa  la  a  channel  of  cons¬ 
tant  cross-section  A  with  ooaponents  of  velocity  0  ~ 
(U  ,0,0)  and  aagnatic  induction  8  “  (0,0,8  )  and 
continuous  electrodes  of  length  a  along  the  x-axis 
on  the  y-walls  at  y  •  0  and  yah.  Thera  is  no  ex¬ 
ternally  applied  electric  potential  and  the 
conductivity  is  scalar  and  uniform  in  a  region 
extending  froa  xaotox-a,  yaOtoy-h  and  z 
a  0  to  z  a  w.  Let  a  be  saallar  than  the  extant  of 
the  giagnatic  field  region  x  along  the  x-axls  so 
that  a  <<  X  and  the  applisd  aagnatic  field  8^  is 
constant  in  ths  region  a. 


To  determine  the  influence  of  Rm  on  the  total 
power  it  is  necessary  to  integrate  Eg.  (6)  over  the 
interaction  region,  which,  because  of  uniformity  of 
the  parameters  involved,  is  egui valent  to  integra¬ 
tion  over  x  in  the  interval  froa  x  •  q  to  x  •  a. 
before  we  attempt  this  integration  note  that  we  can 
obtain  from  the  Maxwell  equations: 

’xB  -  UgJ  (8) 

and  for  ths  purposes  of  this  discussion 

-  IT  ~  “O^y 

whore  the  minus  sign  is  again  due  to  the  geometry. 


At  the  same  time  the  familiar  induction  equa¬ 
tion  can  be  rederived  for  this  special  case  as  fol¬ 
lows:  Froa  Bq.  (9),  by  diffsrentiation  we  obtain 


and,  also  by  diffarantiation, 
wo  obtain 


>x2 

from  Eqs. 


(10) 

(3)  and  (4) 


Currently  favored  expressions  for  the  power 
generated  by  a  siapla,  one-dimensional,  ideal 
generator  are  usually  given  for  the  case  where  the 
magnetic  Reynolds  nuabar  n  is  much  smaller  than 
unity  (Ra  «  X).  These  Xow  Mi  axprosslons  aalca  no 
assumption  concerning  ths  region  of  high  conduct¬ 
ivity  in  the  x-diroction,  a.  In  fact,  they  assume 
that  the  magnetic  field  8  is  a  squars  function 
(i.e.  they  assume  8  X8  ia0;(x£a  and  zero 
eXsawhare) .  Since  for  XeS  Mi,  8^  e  8^,  this  as¬ 
sumption  is  valid.  For  a  continuous  Mm  generator 
the  extent  of  the  high  conductivity  region  can  be 
assumed  to  bs  infinite  (i.e.  o  is  high  in  the 
region  froa  x  ■  0  to  x  ■  •  and  only  the  extant  of 
I  is  limited) . 


ao 

TT 


i) 


(It) 


Since  we  have  assisMd  that  8  and  U  remain 
constant  with  x  for  the  purposes  ^of  this  approx¬ 
imation,  we  obtain 


n 

°K  TT 


o"  ai* 


(12) 


This  partial  differential  equation  has  solutions  of 
the  fora 


.  -Rmx/a 

A  • 


(13) 


In  this  analysis  we  consider  the  extant  of  the 
high  conductivity  region  to  be  limited  to  0  £  x  £  a 
and  a  <  Xg. 


Froa  Ohm's  law  we  obtain,  for  the  given 
gacaatry, 


while  E^'  is  given  by 


E‘  e  E  -  U  B 
y  y  X  2 


(3) 


(4) 


Note  that  the  minus  sign  in  Eq.  (4)  is  a  simple 
consequence  of  the  geometry.  Then  we  obtain 


OB' 


Note  that  Bq.  (13)  yields 


>B 


and 


l!!*  -  -a 

»x*  ^  ■ 


(14) 


(15) 


Mots  that  froa  Eqs.  (12),  (14)  and  (15)  we  obtain 

.0  -  -i-  .  a,  (SS)  (16) 

xxs  x# 

i.e.  Bq.  (13)  satisfies  Eq.  (12). 

Let  the  boundary  conditions  be 


and  the  power  per  unit  voliaa  of  plasma  P  delivered 
to  the  load  (in  the  generator  case)  is  given  by 


♦  ♦ 
p  .  J.E 

where 


J..E 


y  y 


Ky]  f6«> 
(eb) 


The  negative  sign  convention  means  that  power  is 
rsiaoved  from  the  plasma. 


The  magnetic  Reynolds  number  for  the 
geometry  is  given  by 

Rm  i  UgOOj^a. 


given 


(7) 


Bj(0)  s  b”  (17a) 

and 

B^{a)  *  B*  .  (17b) 

Then  observe  tnat  Eq.  (9)  gives 

/*  Jy  dx  .  i  -  -u-^(b;  -  b;)  ('8) 

0 

where  1  is  the  current  delivered  to  the  load  which 
in  terms  of  Bq.  (5)  becomes 
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i  .  -  /*  -  Ky]dx 


(191 


In  tha  liait  in  which  wa  hava  1  <•  O,  tha  aa^^ 
natic  fiald  mat  ba  aquai  Co  tha  appliad  aaqnaeic 


{iald,  .  .  •  Bq-  This  is  aquivalanc  to  daaand* 


ing  that  tha  avaraga  of  tha  adga  fialda  ba  aqual  to 
tha  appliad  aagnacic  fiald,  B^,  tharafora  wa  hava 


(20) 


Anothar  way  of  stating  Eq.  (20)  is  to  say  that 
tha  inducad  aagnatic  fiald  laaadiataly  in  front  of 
tha  plasiaa  Bust  axactly  opposa  tha  corrasponding 
inducad  magnatic  fiald  at  tha  trailing  adga  of  tha 
plassu  or  div  B  ~  O, 


Wa  can  now  procaad  to  solva  for  tha  constant 
•  Fron  Bq.  (13)  wa  obtain  at  x  ~  O 


®z  "  ®0*y  - 


(21a) 


and  at  X  ■  a 


®z  “  *0%  ~  *1*  •  >2ib) 


Tharafora,  froa  Bqa.  (20)  and  (21)  wa  obtain: 


and 


^Vy  -  *1(1  ♦  •  )  -  2Bo 

2Bo(l  -  K  ) 


(22) 


• 


(1  + 


(23) 


Noca  also  froa  Eq,  (18)  wa  obtain 


«z  - 


Rn 

•A  J 

•  u.  -  •  _ ; _ ^ - - -  f(23a) 


j  2Bo(l-Ej(a”"-l) 


0  w 

and  thus  wa  can  obtain  an  axprasslon  for  tha  total 
current  given  by: 


1/w  -  -(2Bg/Mg)  (1-Ky)  .  (23b) 


Finally,  froa  Eqs,  (13)  and  (23)  wa  obtain 


B  «  B  K  +  Rmx/a 

*  °  y  ■ 


24) 


Insarting  Eq.  (24)  in  Eq.  (5)  wa  obtain 


2BQO0^(l-Ky)  ♦  • 


Rmx/a 


d+e’^) 


(25) 


By  dif farantiation  of  Eq.  (24)  wa  obtain 


Hi  .  .  ftn  .  ^Rmx/a  (26) 

3x  * 


'  'j .  >  “j.  • » ■>  ■ 


Tharafora,  insarting  Eqs.  (25)  and  (26)  in  Eq.  (9), 
wa  obtain 


Uo»0x» 


-uoO. 


an  identity •  ( 27 ) 


Tharafora  equations  (24)  and  (25)  also  satisfy  tha 
differential  equation  for  tha  current  density,  Eq. 
(9). 


Wa  are  now  raady  to  integrate  Eqs.  (6)  and/or 
(19)  to  obtain  the  power  as  a  function  of  Bis. 


Sinea  the  electric  fiald  in  the  plasma  due  to 
tha  load  is  constant,  wa  can  writa,  for  the  voltage 


across  tha  load  V^, 


Vg  -  Byh  -  KyUjBoh 


(28) 


Wa  observe  that  the  quantity  [(B^/B^)  -  K.,]  , 


where  only  the  first  term  is  variable,  appears  in 
tha  integrand  of  )»th  Eqs.  (6)  and  (19)  when  they 
are  integrated.  Let  us  therefore  obtain  the 
integral  of  this  quantity.  From  Eq.  (24) 


/*  [(r)  -  I'vl'S*  ■  ^  /* 


0  "0  y 

a 


“  /  [*o*. 

■0  0  '' 


®0  0 
2Bg(l-R..) 


Jhnx/a 


(l+a"“) 


•O'^yJ 


a 

BB'Ba 


/"*  e»”x/a 

0  (l+a“®)  • 


a  ,2Bo(l4E^)a“  .  2Bo(l-K^), 

•  isrrr-  I - - 


d+a"®) 


2Bo« 


R»*B„ 


^RRI  9 


a"®+r 


2«(1-Ky) 

tanh(Rni/2) 


Rm 


(29) 


note  that  this  talres  care  of  the  integration 
w.r.t.  X  only. 


Substituting  Eq.  (29)  in  Che  integral  over 
volume  of  Eq.  (6)  we  obtain  for  the  total  power  to 
Che  load,  P^, 


a,w,h 
J  PdV 
0 


2a(I-iC  ) 

»U?b2  K..  .  - ^ 


'x^O  ‘'y 


UqOUx* 


wh»tanh(Rm/2 ) 


B  2 


-  4  (7^)  IJ,K„(1-K„)  •wh*tanh(Rm/2] 


*  y 


(30) 


The  quantity  wh  appears  when  the  uniform  w.r.t.  y 
and  z  integrand  is  integrated  w.r.t.  y  and  z  and 
evaluated  at  y  •  h  and  z  •  w. 


Wa  can  writa  Eq.  (30)  as  the  power  to  the  load 
par  unit  cross-sectional  area  of  the  generator 
channel,  A  •  wh,  as 
2 


f^/^  -  -  4( 


2u, 


U-K..  (1 


X  y 


Ky)  tanh^Rm/2,  (31) 
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Nota  that  alnca  tanh(m/2)  *  1  aa  Bb  *  •  , 
tha  powar  obtainad  from  a  high  nagnacic  Raynolda 
nuabar  dev  lea  (a.g.  one  with  vary  high 
conductivity)  will  reach  tha  limit 


Eq.  (25),  for  Ra  *  0  yielda 

Jy  -  ‘34) 


P^/A  a  -4(b2/2uo)  (l-KyJ  (32) 

Note  that,  |ince  already  at  Ra  •  3  tanh(RB/2)  > 
0.9,  for  0  •  10  n/sac  and  a  •  1  a,  we  obtain  ^  • 
239  mho  /m.  Tharefore  there  la  an  optimum  conduc¬ 
tivity  for  each  pulsed  MHD  generator  given  by 
O’ 3  S/fi.Ua.  This  expression  can  be  derived  by 
reasoning  that  operation  at  a  higher  magnetic 
Reynolds  number  is  not  required  especially  since 
other  losses  or  costs  can  only  increase  aa  Ibi  is 
Increased.  A  mindless  quest  for  ever  higher 
conductivity  is  not  necessary. 

It  is  now  appropriate  to  note  that  Eqa.  (30) 
and  (31)  both  approach  zero  aa  Dm  «  0 .  In  other 
words,  the  power  to  the  load  vanishes  aa  ito  ■»  0. 
All  other  expressions  are  well-behaved  also.  For 
example  aa  Rm  *  0,  Eq.  (24)  yields 


(33)  and  (34)  are  both  the  low  magnetic 
Reynolds  numbers  approximations  which  yield 


■  (1-K„)a  (35) 


Rn - “0“x‘'y'*  "y- 


o.R™  K„(i-Ky) 


'Up'  X  m  y' 


Note  that  from  Eq.  (36)  we  can  derive  the 
expression  for  power  to  the  load  per  unit  volume  of 
plasma  in  the  generator  given  by 


2Bp(l-K  )-l 

B_  -  B„K..  +  — - i -  (33) 

1  +  1 


’2  “0‘^y 


82  ■  80 


^  .  -O  B^ujKyd-Ky)  (37) 


Equation  (37)  is  the  source  of  some  confusion 
whan  used  for  high  Rm. 


ACCELERATE  I  g  TO  1000  lim/» 
Jy  (A/cm^)  ot  IT  or  Ey  (V/hi)  01 10^  S<1n 


Fig.  1.  Feasible  domain  to  accelerate  1  g  mass  to 
100  )an/s  by  application  of  the  Lorentz 
force  in  an  imposed  magnetic  field  of  1  T 
in  a  linear  MHD  device  of  length  L, 
diameter  D,  and  volume  V.  As  Lorentz 
force  increases,  device  scale  must 
decrease  for  acceptable  aspect  ratios  and 
mass  densities. 


ACCELERATE  100  g  TO  10  km/s 
Jy  iA/cmh  OIIT  or  Ey  (VAn)  ot  lo'^  SAn 


Fig.  2.  Feasible  domain  to  accelerate  100  g  mass 
to  10  )aa/s  by  application  of  the  Lorentz 
force  in  an  imposed  magnetic  field  of  1  T 
in  a  linear  MHD  device. 


9 


«  1984  STD  RESEARCH  CORPORATION 


J 


INTERACTION  PARAMETER  BASED  ON  VELOCITY,  s  J|(jyBj/pu2[  ax 

o 


Pig.  3.  Hap  of  the  Bagnetlc  Reynolda  ntaibar  and  interaction  paraaeter  based  on  velocity  for  some 
significant  HUD  generator  experiments  and  designs  using  as  woricing  fluids  (a)  fossil  fuel 
coaibuation  products,  (b)  solid  propellant  combustion  products,  (c)  explosion  products  and 
noble-gas  driver  gas,  and  (d)  nitrogan  or  argon  arc-jat  plasma. 


INTERACTION  PARAMETER  BASED  ON  PRESSURE,  Ip  3/ j(jyB^)/pjdx 

0 


i 


Pig.  4.  Hap  of  the  magnetic  Reynolds  number  and  interaction  parameter  based  on  pressure  for  the  same 
HHD  generator  experiments  and  designs  given  in  Fig.  1. 
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VERY  HIGH 
INTERACTION: 


INTERACTION: 
STW»g"HgW:iNEAR 
C0UPCIN6  BETWEEN 
FLOW  FIELD,  LOR- 
ENTZ  FORCES  ANO 
JOULE  DISSIPATION 


MODERATE 

INTERACTION: 


asymmetry  ANO 
DISTORTION 


LOW 

INTERACTION: 


GASOYNAMICS 
WITH  SECOND- 
ORDER  CORRECTIONS 


INTERACTION  PARAMETER  TIMES  EFFICIENCY, 

Interaction  levels  and  specific  energy  extraction  of  major  HUD  power  generation  experiments  and 
designs  in  the  U«S.A.  and  the  U.S.S«R.  Solid  points  Indicate  experimental  test  data. 
Interaction  regimes  are  arbitrarily  divided  on  the  basis  of  the  degree  of  MUD  flow  distortion 
in  the  STD  calculations. 


mrcMCTiqi 

(CUSSICM.  Msorwmics  mni  ism  coMsaigiiu 


NOPOMIT  IHTOIMTtOi 

(stsiirtcMT  nm  astnstit  m  distmtisi  ns  to  loaoin  fodccs) 


ftmwmmtmmW} 


tmt  IPS!  STSTK  TWMWUnSW  vm  c 


Development  of  velocity,  temperature,  and  currant  under  conditions  of  AEDC/HPOE  Run  MI-006-0 U. 
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static  tonperature  and  transverse  current  distributions  in  the  ABDC/HPDE  with  hot  electrode  walls  (Tw,max  *  1700  K)  and 
peak  magnetic  induction  in  the  channel  of  (a)  3  T,  (b>  4  T#  and  Ic)  S  T* 
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Pig.  11.  Total  power  density  relative  to  the  core 
power  density  st  three  stations  In  a  1000 
MW  (theraal  Input)  coal-fired  MHD 
generator  with  circular  cross  section,  as 
a  function  of  the  angle  subtended  by  the 
electrodes.  Also  shown  (dashed)  arc  the 
Idealised,  constant  property  results  of 
F.J.  Flshaan,  J.  Advan.  Energy  Conv., 
4:223-236  (1964). 


ANODE 


Fig.  13.  Current  Streamlines  In  the  cross  section 
at  the  S  a  station  In  Fig.  11  for  an 
electrode  subtended  angle  of  ^0°..  Peek 
trsnsverse  current  density  0.3  A/ca^. 


Fig.  12.  Current  streamlines  In  the  cross  section 
for  representative  cases  from  Fig.  11. 


ANODE 


Fig.  14.  Streaallnes  of  the  secondary  flow  velocity 
corresponding  to  Fig.  13.  Peak  transverse 
velocity  coaponents  are  9Z  of  the  aean 
axial  velocity. 
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